Universal potential barrier penetration by initially confined wavepackets 
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The dynamics of an initially sharp-boundary wavepacket in the presence of an arbitrary potential 
barrier are investigated. It is shown that the penetration through the barrier is universal in the 
sense that it depends only on the values of the wavefunction and its derivatives at the boundary. 
The dependence on the derivatives vanishes at long distances from the barrier, where the dynamics 
are governed solely by the initial value of the wavefunction at the boundary. 

PACS numbers: 03.65.-w, 03.65.Nk, 03.65.Xp. 



Introduction The tunneling phenomenon is a funda- 
mental problem in quantum mechanics. Traditionally, it 
is solved and presented in most textbooks for a particle 
with a given energy and momentum (see, for example 
0). Although this treatment is usually sufficient to 
demonstrate the main features of tunneling, it raises a 
serious difficulty, for the uncertainty principle suggests 
that the particle is instantaneously distributed over 
the entire space. Clearly, this analysis implies that the 
particle was and always will be at both sides of the 
barrier, and therefore the classical meaning of tunneling 
is not clear. This problem is one of the reasons that 
researchers investigated the dynamics of wavepackets, 
usually Gaussians, in the process of tunneling 0. 
Although this treatment seems more intuitive, whenever 
the initial wavepacket is analytic there is always initially 
a tail at the other side of the barrier. To prevent this 
problem, the initial wavepacket must be singular and to 
have a non-zero value only at a finite region in space. 
Clearly, tunneling is just one specific case in a much 
broader family of scattering problems, and all that was 
said above can be generalized to transmission through 
any arbitrary potential. 

During the last few years, due to the recent development 
in femto-second pulses lasers j3lQ, optical tweezers and 
atom cooling and trapping Uly, an interest in Moshin- 
sky's shutter problem]]} was rejuvenated pi. I9L IToL fTTl Il2| . 
Moshinsky was the first to investigate the dynamics of 
an initially singular wavefunction. The singularity of the 
initial wavefunction modeled a fast shutter. The ability 
to localize particles (usually atoms) by laser beams 
and then to release them instantaneously increases the 
feasibility to simulate the shutter flallij . 
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In this paper we derive a generic formula for the 
Schrodinger dynamics of an initially singular wavefunc- 
tion (which vanishes at half space) in the presence of a 
finite-width potential. 

It should be stressed that even if the initial wavefunc- 
tion is a smooth function (rather than singular) with a 
transition length scale of e, then the results presented 
here are valid provided the measurements are done at 
distances shorter than tH/2me (see ref . [TT| ) . 
The general case A wavefunction is confined to one 
side of a potential barrier. The barrier is located around 
L > and vanishes for x < and x — > 00. The initial 
wavepacket, on the other hand, vanishes beyond x > 
(see Fig.l). Therefore, initially, there is no overlap be- 
tween the wavepacket and the barrier. The dynamics are 
governed by the Schrodinger equation 



dx 2 



i> + V (x - L) ip = i 



dip 



(1) 



where V (x) is the potential barrier, 
units H = 1 and 2m = 1 are used. 



Hereinafter, the 



Since the initial wavefunction has a singularity at x = 
it can be written as[ll|, [l2j ip (x, t = 0) = / (x) 8 (—x), 
and the solution at t > and x > L is then 



ijj (x,t) 



dkda F{q) _JT^)_ e ^ { . kx _. kH) 



2tt 



k - q + iO 



(2) 



In this equation F (q) = (2n) 1 J dxf (x) exp (— iqx) is 
the Fourier transform of / (x) and T (fc) is the transmis- 
sion coefficient of the barrier for plane- wave with momen- 
tum k. 

The properties of the barrier and the initial wavefunction 
can be separated by defining 



1 f iT (k) 1 9 \ 

ld (q x, t) = —— / dk- exp (ikx — ik t) 

2n J k — q + iQ 



(3) 
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FIG. 1: scattering of an initially singular wavepacket over a FIG. 2: An illustration of the four scenarios under study. A) 
barrier. In the figure only the real part of the wavefunction is delta function potential B) Arbitrary smooth opaque barrier 
plotted. C) Rectangular barrier and D) Reflectionless potential well 



and hence 



ip(x,t) = / dq<p(q,x,t)F(q), 



(4) 



which, by expanding ip (g, x, t) in orders of q, Eq. @ 
can be rewritten in terms of f(x). We thus obtain the 
generic solution of the problem: 



ip(x,t) = ip 



-t^,s,t)/(i> 



This generic expression should be interpreted as the in- 
finite expansion 



(-i) n d n d n 



ip(0,x,t)f(0)-i<p'(0,x,t)f(0) + ---, (6) 

where <p' (0,x,t) = dp (q, x, t) /dq\ q=0 and /' (0) = 
d f (x) /dx\ x=0 . 

Eqs. © and I© are the main result of this paper. It 
is evident from Eq.© that the wavefunction dynam- 
ics depend only on the initial wavepacket values at the 
boundaryir = 0. It is a universal formula in the sense 
that any wavefunction / (x) and any potential can be 
substituted in it to find the exact dynamics beyond the 
potential at any t > 0. In particular, when / (x) = C is 
a constant, then the solution is simply 

fl>(x,t)=C<p(0,x,t). (7) 
Eq.© can also be written as the following expansion 



^(x,i) = E a «/ (n) (°)> 



(8) 



71=0 



where /(") (0) = d n f (£) /<9f l | ?=0 are the initial func- 
tion's derivatives and the coefficients are 



27m! 



dk- 



T(k) 



(k + iO) 



j exp [ikx — ik 2 i) . (9) 



These coefficients can be written also as an expansion 
by expanding T (k) j [k + i0) n+1 around the momentum 
x/2t (the dominant term in the stationary phase approx- 
imation in long times and arbitrary^), i.e., 



T(fc) _ 



(k-x/2t) m d m [T(q)/q n+1 ] 



(10) 



q=x/2t 



(5) By substituting Eq. (jTU|) in Eq.© and Eq.© one obtains 



V,( a;) t) = -e X p(^)^6 n /W(0), (11) 



n=0 



where the coefficients arc 



n\ ' (2to)! 

m=0 v ' 



r( m + r)(ity {2m+1)/ \ (12) 



and 



Sn,2r. 



_ d A 



>)k- 



[T (k)k- 



lnt . This result also 

J k—x/2t 



suggests that at large distances x jt — > oo 

*(».•) ~7-y^— /(0), (13) 

and the particles density decays like p = |?/;(£,t)| 2 ~ 
\T{x/2tyf(0)\tx- 2 /Tr, which agrees with the free 
resultsjn| for T = 1. Moreover, for every finite potential 
T (k — > oo) — > 1, which means that at very short times 
(with respect to the distance, i.e., t/x 2 « 1, although 
it was derived for large t) Eq. Ijl3|l is reduced to the free 
propagation for every potential, which is consistent with 

Specific examples and applications In this section we 
will demonstrate these derivations in a few scenarios for 
different potentials. 

The delta function potential - an exact solution. 

One of the simplest scenarios of a potential barrier pene- 
tration is the delta- function potential (see Fig.2A). This 
problem was discussed in detail for an initially singular 



3 



wave function in [12J , and is presented here only to illus- 
trate the validity of the above derivations. In this case 
the Schrodinger equation reads: —§^tp + X5 (x — L) tp = 

I 



i^f- The transmission coefficient in this case is well 
known T (k) = k/(k — i\/2), in which case by Eqs.© 
and 1^1 the solution beyond the barrier is 



ip{x> L,t) 



2 /At 



d/d^ + A/2 



, x d 



r 



A/2 



d/d£ + A/2 



,i y .A 
tt\ - — i— 

2t 2 



f(0\ 



(14) 



This expression obviously agrees with [12j , and illustrates 
the fact that the dynamics are fully described by the 
values of the initial wavefunction at the singular point. 
At short times, i.e., t « x 2 , it can be expanded (see 
0) 



fit e" 2 / 4t 

7T X 



(15) 



-t)/(o) 

x 



2t 



iftfW (0) + O(x~ 2 ) 



This expression is, of course, consistent with Eqs. © 
and 10. 

A special case occurs when A = 2/W (0) // (0). In this 
case the leading terms have no dependence on the barrier: 



ip (x,t) 



it e " 2 /4t 



7T X 



[/(0) + O(x- 2 )]. (16) 



That is, in the leading terms, the influence of the bar- 
rier was eliminated by changing the derivative of the 
initial wavefunction only at the singularity point for 
\ftjx « 1. In general, each one of the terms in the 
expansion Ijl5(l can be eliminated (except for the first 
one) by tailoring the derivatives of the wavefunction at 
the singular point. 

An arbitrary smooth opaque barrier. In the case of 
an opaque barrier with smooth edges (see Fig. 2 B) one 
can use the WKB approximation to evaluate the trans- 
mission coefficient of the barrier. 

If the potential function is V (x) then [l| T (fe) = 

(29+1/20) . where 9 ( fc ) = expC/o^t^)^ 2 ] 172 )- 
Therefore, from Eci. (fTT|) 



itf(0) 2exp(ix 2 /4t) 

7T x [26(x/2t) + l/29(x/2t)}' 



(17) 



Eq. i|17|) is valid, of course, only when the turning points 

exist, i.e., when (x/2t) 2 < max{V" {x)}. 

For very opaque barrier, i.e., di][V (77) — (x/2t) 2 ] 1 / 2 ^> 

1, 

itf(0) e lx2 / 4t Q fx\ 



^^~^^— d {2t) 

The rectangular potential barrier. In the case of the 
rectangular potential barrier (see Fig.2C), with potential 



height and length Vq and L respectively, the transmission 
coefficient is[l[ 



T(k) 



exp (—2ikL) 



cosh (2kL) + i (e/2) sinh (2kL) 



(18) 



, where n = VVo — k 2 and e = njk — kj ' k. 

At short times, the particles that penetrate the barrier 
are extremely energetic and therefore they pass the bar- 
rier almost unaffected by it. 

In this case k = ix/2t, e = i 1 + 2 (Vot 2 /x 2 ^ 2 and 



T = 1 - 2i (V t 2 /x 2 ) 2 sm (xL/t) exp (ixhjt) 
which yields by Ea. (|13|l 

e ix 2 /4tf (0) ^/ttj^/x 



tp (x,t) 



1 - 2i {V t 2 /x 2 ) sin (xL/t) e ixL / 1 



,(19) 



with weak resonances at xL/t = rmr, (where m is an 
integer). 

A Reflectionless potential. One of the peculiar poten- 
tial examples is the reflectionless potential^, 0|. For 
this potential (see Fig. 2D) the absolute value of the 
transmission coefficient is always |T(A;)| = 1. Therefore, 
the reflection coefficient is zero for every incoming plane 
wave; however, the transmission coefficient suffers from 
dispersion, which deforms the initial wavepacket. 
It is well known that the potential V (x) = 
— 2a 2 sech 2 (ax), whose width goes like ~ a~ x and its 
depth like a 2 , belongs to a family of reflectionless po- 
tentials [H El E3- 

This potential, as opposed to the previous three, is not 
completely localized in space. That is, the initial wave- 
function "feels" the potential at any distance. Therefore, 
one should place it at a large distance from the edge of 
the initial wavefunction. 
Its Schrodinger equation is then 



dx 2 



ip — 2et 2 sech 2 [a (x — L)] ip = i 



~dt' 



(20) 



where L >> a 1 is the new location of the potential well. 
This equation has an exact reflectionless solution (which 
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corresponds to the incoming plane wave ip (x — > — oo) ~ 
exp (ikx) + R (k) exp (— ikx) ) 

ip (x) = [k/a + i tanh (ax)] (k/a — i) 1 exp (ikx) . (21) 

For x — -> oo Ea. (|21[l can be written ip (x) ~ 
T (k) exp (i&x) with the simple transmission coefficient 
T (k) = |±g. According to Eqs.©, 

if(q,x>L,t) = ± J dk k _l +i0 k +_T a exp (ikx - ikH) , 



which can easily be written as [12j 

ip(q,X,t) = ip free (q,X,t) - (22) 
q [V/ree (<?, t) - <p/ ree (a, X, t)] , 

where (ff r ee (q, x, t) is the free (without the barrier) solu- 
tion (pf ree (q, x, t) = |w U/it (§: — q)] exp (ix 2 /2t) . The 
solution according to Eq.JSJ) is finally 



ip(x,t) 



free [ i 7 ^7 1 



2ai 



d/di + a 



free 

r 



/(Ol 



(23) 



In a particular case, where there is a discontinuity only 
in the function, and not in its derivatives, i.e., if the initial 
wavefunction looks like ip (x, t = 0) = / (0) 9 (—x), then 



ip (x, t) 




^/(O)expfe 



2i 



(24) 




V V2i 



Again, at short times t << x/a the free solution is 
retrieved (the two expressions on the right cancel each 
other). 

Conclusions and summary 

We have presented a generic formula, which de- 
scribes the solution of the temporal propagation 
through an arbitrary potential of an initially sin- 
gular wavefunction. It was shown that when the 



initial wavefunction vanishes at the entire half space, 
i.e., ip (x, t = 0) = / (x) 6 (— x) the wavefunction for 

, where 



ik 2 t) is a 



every t is ip (x,t) = ip (-i^, ar, *) / (0 
f(q,x,t) ee jL / dk k _g +M T (k) exp (ikx - 

function, which depends only on the transmission 
coefficient of the barrier T (k). This expression was also 
generalized to an initially compact support wavefunction. 



Particularly, it is shown that at very large distances 
the solution can be approximated by the generic expres- 
sion for any initial function ip (x,t = 0) = f(x)8(—x) 



the general solution is ip (x, t) ~ T (||) J ^ 



/(0). 



The authors are indebted to Miriam Schler for her net- 
working efforts. 



[1] 15. E. Merzbacher, Quantum Mechanics, (John Wiley 

and Sons, New- York 1970) 
[2] A. P. Stamp and G.C. Mcintosh, Am. J. Phys. 64, 264 

(1996); L.I. Schiff, Quantum Mechanics, 3 rd edition 

(McGraw-Hill Publishing, Singapore, 1986). 
[3] G.G. Paulus, F.Lindner, H. Walther, A. Baltuska, E. 

Goulielmakis, M. Lezius, and F. Krausz, Phys. Rev. Lett 

91, 253004 (2003) 
[4] M. Hentschel, R. Kienberger, Ch. Spielmann, G.A. Rei- 

der, N. Millosevic, T. Brabec, P. Corkum, U. Heinzmann, 

M. Drescher, and F. Krauz, Nature 414, 509 (2001) 
[5] C.E. Wieman, D.E. Pritchard, D.J. Wineland, Rev. Mod. 

Phys., 71, S253 (1999) and references therein. 
[6] K.C. Neuman and S. M. Block, Rev. Sci. Instrum. 75, 

2787 (2004) and references therein. 
[7] M. Moshinsky, Phys. Rev. 88, p. 625 (1952) 
[8] S. Brouard and J.G. Muga, Phys. Rev A, 54, 3055 (1996) 
[9] S. Godoy, Phys. Rev A 65 042111 (2002) 
[10] F. Delgado, J.G. Muga, and A. Ruschhaupt, Phys. Rev 

A 69, 022106 (2004) 



[11] E. Granot and A. Marchewka, Europhys. Lett. 72, 341 

(2005) ; A. Marchewka and Z. Schuss, Phys. Rev A 61, 
052107 (2000) 

[12] E. Granot and A. Marchewka, Phys. Rev A, 73, 032111 

(2006) ; A. Marchewka, E. Granot and Z. Schuss, ac- 
cepted for publication in Optics and Spectroscopy 

[13] 13. P. Szriftgiser, D. Guery-Odelin, M. Arndt, and J. 
Dalibard, Phys. Rev. Lett. 77, 4 (1996) 

[14] 14. C. Fort, L. Fallani, V. Guarrera, J. E. Lye, M. Mod- 
ugno, D.S. Wiersma, and M. Inguscio, Phys. Rev. Lett. 
95 170410 (2005). 

[15] 16. L.D. Landau and E.M. Lifshitz, Quantum Mechan- 
ics, 3 rd ed. Pergamon Press, Oxford 1985, sec. 25 and 
problems therein. 

[16] 17. N. Kiriushcheva and S. Kuzmin, Am. J. Phys., 66, 
867 (1998) 

[17] 18. B.Y. Tong and N. Kiriushcheva, Phys. Lett. A 229 
49 (1997) 



